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In this talk we present the extension of the antenna subtraction method to include initial states 
containing one hadron at NNLO. We sketch the requirements for the different necessary subtrac- 
tion terms, and we explain how the antenna functions are integrated over the appropriate phase 
space by reducing the integrals to a small set of master integrals. Where applicable, our results 
for the integrated antennae were cross-checked against the known NNLO coefficient functions for 
deep inelastic scattering processes. 
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1. Introduction 

Final states containing hadronic jets are produced at large rates at high energy particle collid- 
ers. Owing to their large production cross sections, various jet observables can be measured to a 
high statistical accuracy. However, experimental data on these observables are often so precise that 
meaningftil precision studies must rely on theoretical predictions that in perturbative QCD requires 
corrections at next-to-next-to-leading order (NNLO). 

At NLO and NNLO the contribution from real and virtual corrections are separately divergent 
and, while infrared singularities from purely virtual corrections are obtained immediately after inte- 
gration over the loop momenta, their extraction from purely real or mixed real-virtual contributions 
is much more involved. In the latter case the singularities become expUcit only after integrating the 
matrix element over the appropriate phase space. Since the integration is in most of the cases not 
feasible analytically, one need to extract infrared divergencies at the integrand level and construct 
subtraction terms which should satisfy the following two conditions: (A) they should approximate 
the real radiation matrix element in all singular limits, and (B) they should be sufficiently simple 
to be integrated analytically over a section of the phase space that encompasses all regions cor- 
responding to singular configurations. In the contribution to an m-jet cross section at NLO and 
NNLO the subtraction terms are added and subtracted in the following way: 

do-NLO = / (dc7i^o-do-|Lo)+ / dc7^L0+ / dcT^LO (1-1) 

dClNNLO = / (d<7NNLO - da^NLo) + / ( d^^NNLO " ^a^NLO ) 

/ d(7^o+/ ^<^o+l daj^j^io. (1-2) 

where Aaf is the real radiation contribution, Aaf and da^^^Q the subtraction terms of the real 
and real- virtual contributions respectively and da^'" is the virtual contribution at n loop (/ =NLO, 
NNLO). 

In the last decades several subtraction methods were developed at NLO [1-4] and NNLO [5- 
10]. One of these methods is the so-called antenna subtraction, which was derived at NNLO in [1 1] 
for partons only in the final state. The anteima subtraction formafism constructs the subtraction 

terms from antenna functions. Each antenna function encapsulates all singular limits due to the 
emission of one or two unresolved partons between two colour-connected hard radiator partons. 
This construction exploits the universal factorization of matrix elements and phase space in all unre- 
solved Umits. The anteima functions are derived systematically from physical matrix elements [12]. 
For processes with initial-state partons, anteima subtraction has been fully worked out only to NLO 
so far [13]. In this case, one encounters two new types of antenna functions, initial-final antenna 
function with one radiator parton in the initial state, and initial-initial antenna functions with both 
radiator partons in the initial state. In this talk we present very recent results on the derivation of all 
NNLO initial-final antenna functions and their integration over the appropriately factorized phase 
space. Finally we show how our results could be cross-checked for many cases against the known 
NNLO coefficient functions for deep inelastic scattering. The initial-final antenna functions form 
part of the full set of anteima functions needed for NNLO calculations of hadron collider processes, 
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and are, together with the already known final-final antenna subtraction terms, sufficient for NNLO 
calculations of jet observables in deeply inelastic lepton-hadron scattering. 

2. Initial final antenna subtraction at NNLO 

In the case of one paiton in the initial state the subtraction term at NLO has the following form 

m+l "^"'+1 

km+l\Xp,r)^ jIi ''(ki, . . . ,K[[, ■ ■ ■ ,^m+l) ) (2-1) 

j 

where / labels the hard radiator with momentum p in the initial state. The additional momentum r 
stands for the momentum of the second incoming particle which can either be colored or colorless. 
This contribution has to be appropriately convoluted with the parton distribution function fi. The 
tree antenna X^.j^, depending only on the original momenta p, kj and k^, contains all the configu- 
rations in which parton j becomes unresolved. The ni-parton amplitude depends only on redefined 



on-shell momenta k^,.. . ,Kk, ■ ■ ■ , and on the momentum fraction x. The jet function, jin'\ in ( |2. l] ) 
depends on the momenta kj and kj^ only through K^. 



As eq. ( |1.2D shows, at NNLO two types of contributions to m-jet observables require sub- 
traction: the tree-level m + 2 parton matrix elements (where one or two partons can become unre- 
solved), and the one-loop m + l parton matrix elements (where one parton can become unresolved). 
In da^^^Q, we have to distinguish four different types of unresolved configurations: 

(a) One unresolved parton but the experimental observable selects only ni jets; 

(b) Two colour-connected unresolved partons (colour-connected); 

(c) Two unresolved partons that are not colour connected but share a common radiator (almost 
colour-unconnected) ; 

(d) Two unresolved partons that are well separated from each other in the colour chain (colour- 
unconnected). 

Among those, configuration (a) is properly accounted for by a single tree-level three-parton an- 
tenna function like used already at NLO. Configuration (b) requires a tree-level four-parton an- 
tenna function (two unresolved partons emitted between a pair of hard partons), while (c) and (d) 

are accounted for by products of two tree-level three-parton antenna functions. On the other side, 

vs 1 

the one-loop single unresolved subtraction term dOj^j^^o ^lust account for three types of singular 
contributions: 

(a) Explicit infrared poles of the virtual one-loop (m + l) parton matrix element. 

(b) Single unresolved limits of the virtual one-loop (m+l) parton matrix element. 

(c) Terms common to both above contributions, which are oversubtracted. 
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For all these cases the detailed form of the subtraction term is given in [14]. The only genuinely 
new ingredient appearing at NNLO are the four-parton initial-final antenna function ^fjj^i, which 
can be obtained by crossing the corresponding final-final antenna functions, and one-loop three- 
parton initial-final antenna function Xfjp. which can be obtained by crossing from their final-final 
counterparts, listed in [11]. Both have to be integrated over the appropriate phase space. 

3. Integration of initial-final antenna functions at NNLO 

The initial-final antenna functions all have the scattering kinematics q + pi ^ Pi + Pii+Ps), 
where 

q^ = -Q^<0, pj = 0, Z=x^, p] = pI=pI=0, 

2q-pi 

and p-i is present only for the NNLO real radiation antenna functions. Thus, integration over the 
final-state two-parton or three-parton phase space yields a result which depends only on and z- 
The NNLO double real radiation antenna functions ^,-'^^/ have to be integrated over the in- 
clusive three-parton final state phase space. The NNLO one-loop single real radiation antenna 
functions X/y^ are integrated over the inclusive two-parton final state phase space, and over the 
loop momentum. For both types of integration, we first expressed all phase space integrals as loop 
integrals with cut propagators [15], and then employed the by-now standard technique of reduction 
to master integrals using integration-by-parts (IBP, [16]) and Lorentz invariance (LI, [17]) identi- 
ties among the integrals of any given topology. After carrying out the reduction, we found nine 
master integrals for the NNLO double real radiation antenna functions and six master integrals for 
the NNLO one-loop single real radiation antennae. All the masters integrals are shown in Figure |l|. 
The convention for naming the master integrals follows the labelling of the numerators, i.e. 

,^,j^,^^ j\^P}mm5\q + p,-p,-p,-p,), where [d;,] = (;,^) . (3.1) 

For the explicit definition of the propagators D, used in Figure [l] we refer to [14]. 

All masters, except /[1, 2, 4, 5], have been computed by direct integration and by the differential 
equations method, supplemented, where necessary, by a direct calculation at x = 1 after factorizing 
the leading singularity. Where appropriate, we compared our results to the expressions in the 
appendix of [18], finding full agreement. The explicit results of the master integrals up to the 
needed order in e as well as the integrated initial-final antenna functions can be found in [14]. 

4. Rederivation of NNLO coefficient functions 

Being derived from physical matrix elements, the integrated antenna functions can be com- 
pared to results from literature for inclusive cross sections or coefficient functions, as was done 
previously for the final-final antennae in [12,22]. In the case of the initial-final antennae, we can 
compare the integrated quark-antiquark antennae and gluon-gluon antennae against NNLO correc- 
tions to deep inelastic coefficient functions known in the literature. The former ones can be checked 
against DIS structure function calculations [18] whereas the latter can be compared to the (/)-DIS 
structure functions computed in [19,20]. These structure functions are obtained in an effective the- 
ory with a scalar coupled to the square to the gluon field strength tensor. For example, the leading 
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I[l,3,4,6] I|2.3.5,6] I[l,2,4,5] 



Figure 1: Master integrals for the phase space integration of the double real tree level initial-final antennae 
at NNLO (left), and for the loop plus phase space integration of the one loop initial-final antennae at NNLO 
(right). The double line in the external states represents the off-shell momentum, q with = —Q^, the single 
one is the incoming parton. All internal lines are massless. The cut propagators are the ones intersected by 
the dotted line. 



colour piece of the two-loop gluon initiated structure function can be written as the following Unear 
combination of antennae: 



(2) 



If- 



If 

l,R 



where ^g^ggg is the integrated tree level gluon-gluon double real radiation antenna, ^g',gg the in- 
tegrated one-loop gluon-gluon antenna and F^^^ and F^^^ are respectively the one- and two-loop 
coefficients of the gluon form factor given in [21]. An explicit expression for the two-loop quark- 
and gluon-initiated structure functions can be found in [14,20]. The explicit linear combinations of 
antennae reproducing the different color contributions of the coefficient functions are given in [14]. 

The quark-gluon antennae, derived from neutralino decay, cannot be associated to any physical 
process and only the deepest pole structure could be checked against a combination of Altarelli- 
Parisi spUtting functions. 



5. Conclusions 



In this talk, we presented the extension of the NNLO antenna subtraction formaUsm [11] to 
include initial-final antenna configurations, where one of the hard radiator partons is in the ini- 
tial state. Furthermore a highly non-trivial check of our results was performed by rederiving the 
two-loop coefficient functions for deep inelastic scattering. The subtraction terms presented here 
allow the construction of a parton-level event generator program for the calculation of NNLO cor- 
rections to jet production observables in deeply inelastic electron-proton scattering. Moreover, the 
initial-final antenna functions derived here are an important ingredient to the calculation of NNLO 
corrections to jet observables at hadron colUders, which will be possible once the computation of 
the initial-initial antenna configurations will be accomplished [23]. 
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